STRUCTURE OF CORONAL NEUTRAL SHEETS G, W, Pneuman 


A qualitative model for the structure of the neutral sheet lying along the axis of coronal ABSTRACT 
streamers is developed. The internal topology of the sheet is that of extremely thin 
magnetic tongues greatly distended outward by the solar wind expansion inside the sheet. 

Due to finite conductivity effects, expansion is taking place across the field lines but is 
retarded relative to the external flow by the reverse jXB force. The sheet thickness is 
determined by three considerations: the electrical conductivity ’ that specifies the magni- 
tude of the gradients in field strength, the expansion velocity that stretches the field lines 
outward decreasing the sheet thickness, and finally, the lateral pressure balance that limits 
the approach of the oppositely directed external field toward the neutral plane. 

If o is the electrical conductivity, the sheet thickness is shown to be proportional to 
a~l/3. For an electron conductivity evaluated perpendicular to the internal field in the 
sheet, the thickness is of the order of 100 km in the inner corona and 10,000 km at 1 
AU. Microturbulence and instabilities are expected to yield dimensions greater than these 
theoretical values since these effects tend to reduce the “effective” conductivity. 


INTRODUCTION 

It is clear that a great number of neutral sheets should 
exist in the solar corona, considering the numerous 
polarity reversals in the observed photospheric magnetic 
field. The intersection of these sheets with the ecliptic at 
1 AU can be observed with satellites and are normally 
referred to as “sector boundaries.” However, the detec- 
tion of structures as thin as these in the inner corona 
would be difficult through normal photographic tech- 
niques. Since an eclipse photograph of the inner corona 
reflects an integration of electron density over the line 
of sight, these sheets, even though denser than their sur- 
roundings, would be essentially transparent unless they 
were not only planar but also had their plane coinciding 
excactly with the line of sight direction. 

There seems to be an increasing amount of indirect 
observational evidence that such sheets may be the 
tracks along which high energy particles released during 
flares travel from the sun to the earth. For example, the 
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most favored theory of type III radio bursts proposes 
that this emission is produced by a cloud of outward 
traveling electrons exciting plasma oscillations at succes- 
sively higher levels in the corona. It has been proposed 
[Wild, 1964] that these electrons move out along the 
neutral sheets situated along coronal streamers. In 2 iddi- 
tion, Bumba and Obridko [1969] have suggested that 
proton flare activity associated with Bartels’ active longi- 
tudes occurs in the neighborhood of the sector bound- 
aries of the interplanetary magnetic field. 

In this paper, we describe a physical model of the 
expected structure of such coronal neutral sheets. 
Though the approximations are crude, it is hoped that 
the model will provide physical insight as to the grossest 
features that would be expected to be associated with 
these structures. 

DESCRIPTION OF THE MODEL 

Since the structure of a neutral sheet essentially is pro- 
duced by finite conductivity effects, the simplest way to 
visualize it is to first consider a coronal streamer for the 
case a = oo, introduce a small amount of resistivity to the 
plasma, and examine the consequences. A typical helmet 
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streamer (fig. 1) consists of a region of closed mag- 
netic loops in the corona with open field lines adjacent to 
and above the loops, and is commonly observed during 
solar eclipses. For a = a neutral point exists at the top 
of the closed loops, and true sheet currents of zero 



insert shows a blown up section of the neutral sheet 
along the steamer axis. 

dimension exist both above the neutral point and below 
it between the open and closed regions [Pneuman and 
Kopp, 1971]. These lower sheet currents exist because 
expansion occurs along the open field lines but not along 
the closed ones, resulting in a pressure discontinuity 
between the two regions that must be balanced by a 
jump in magnetic field strength. The location of the 
neutral point is determined by the condition that the 
pressure differential between the closed and open re- 
gions at that point is exactly balanced by the magnetic 
pressure in the open region; that is. 

Bit 

Pd-Psru 

where cl refers to the closed region and st to the 
streaming region. 

We now investigate the consequences of introducing a 
small resistivity into this model. Obviously , the first will 
be reconnection of open field lines at the neutral point. 
The outward part of the reconnected field will just be 
expelled outward along the sheet, whereas the inward 
part will become a new closed loop attached to the solar 
surface. Consequently, a hydrostatic pressure distribu- 
tion will be set up along this loop. However, the new 
loop will lie outside that for which equation (1) is satis- 
fied, and the pressure differential will be larger than the 
external magnetic pressure. As a result, the loop will 
expand outward along the sheet and, were the conduc- 
tivity infinite, it would continue expanding to infinity 
(since the field and material are frozen). Hence, this loop 
would eventually become a new open field line and the 
configuration would revert to that which existed before 
reconnection took place. If, on the other hand, o is 


finite, a small amount of relative motion between the 
field and gas is permitted and the loop will not be ex- 
pelled to infinity but will only be pulled outward to 
some location where the diffusion of field backward 
through the material is just balanced by the^ outward 
convection of the field lines. When a steady state is 
achieved, the neutral sheet will have a finite thickness 
and will contain loops of weak magnetic field attached 
to the solar surface with the plasma expanding outward 
across the loops. The resulting configuration might look 
something like that shown at the top of figure 1. Since 
the relative diffusion of the field through the gas is ex- 
pected to be very small for coronal conditions, these 
internal loops would be extremely distended and, con- 
sequently, the transverse field there should be very 
small. 

If this physical picture of the neutral sheet is 
adopted, the sheet thickness will depend on three 
factors: 

1. The magnitude of the electrical conductivity. This 
determines how large the gradients in magnetic field 
strength can be within the sheet. As one would ex- 
pect, the sheet dimension should vary inversely with 
the conductivity. 

2. The expansion velocity within the sheet. Large veloci- 
ties will stretch the field lines more and decrease the 
sheet thickness. 

3. The lateral pressure balance between the sheet and its 
surroundings. This consideration will limit the ap- 
proach of the external field towards the neutral 
plane. 

In the next section, we employ these concepts to 
develop a first-order model of coronal neutral sheets. We 
employ the MHD approach throughout. Since these 
structures are expected to be thin, one must ascertain 
whether collisions will be frequent enough to make this 
approximation valid. However, since the sheet dimen- 
sions are not known a priori, \i seems most fruitful to 
investigate this question in the light of the properties of 
the resulting solution. 

MAGNETIC FIELD CONFIGURATION 

Consider a coronal streamer, situated at the solar 
equator, which is axisymmetric with respect to the rota- 
tional axis. If Q is the polar angle; then, the neutral sheet 
will lie along the plane 6 = tt/2. As a first approximation, 
the sheet is assumed radial in cross section and the fluid 
velocity inside is also taken to be radial.^ In lirie with 
the arguments of the previous section, the magnetic field 


^Throughout this work we omit the effects of solar rotation. 


56 



will have both radial and transverse components. Hence, 
the velocity and magnetic field vectors are given 
respectively by 

V = and B ~ 

where er and bq are unit vectors in the radial and polar 
directions. 

For steady state conditions, the expression for con- 
servation of magnetic flux and the induction equation 
are 


7 • B = 0 



(2) 

(3) 


where a is the electrical conductivity . In deriving equa- 
tion (3), we have taken a to be a scalar. This would not 
be the case in the presence of a strong magnetic field, 
but the field inside the sheet is expected to be weak and 
the reduction in conductivity due to the field may not 
be great. This question will be considered again later so 
that actual conditions may be more realistically assessed. 
We merely consider here that the conductivity is a func- 
tion of temperature only; we also assume isothermal 
conditions such that a is a constant. 

Subject to the above conditions, equations (2) and (3) 
can be solved in the limit » 1 to yield expressions 
for both the radial and transverse components of the 
magnetic field. We merely tabulate the results here. 


Br 
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P„(COS e^y) 


(5) 

Here,x -r/ro, where Tq is the radial distance to the top 
of the helmet in the streamer. This is where the sheet is 
assumed to originate and is where the neutral point 
would be located if the conductivity were infinite. is 
the magnetic field strength at Tq just outside the sheet: 
u = Vr/V^ (Ys being the sound speed); Rem is a magnetic 
Reynolds number given hy = AkotoVs', Pn{co%B) is 


a Legendre polynomial of order n \ and 0^ is the angle 
representing the sheet boundary. The prime denotes 
differentiation with respect to 0, (Note that for a 
coronal temperature of 1.0X10^ °K and 

~ 4.5X10^"^, This is an enormous number and 
shows just how effectively the field and material are 
frozen in the corona.) To satisfy the two conditions that 
By^O for 0 = 7 t/ 2 and Bq = 0 for 0 = 0^, we require that 
n he odd and also a root of the equation 

P^(cos0^,) = O (6) 

Since n characterizes the variation of quantities across 
the sheet, it will be a large number inversely propor- 
tional to the sheet thickness. For example, if n » l ,we 
can write 

” R" * j) « * I ] 


Condition (6) then reduces to 



_ 

n-2djj 


and since nl2 ~ 0j - d/r where 6 is the sheet thickness, 
we have 



(7) 


CURRENT DENSITY, PRESSURE DISTRIBUTION, 
EXPANSION VELOCITY, AND SHEET THICKNESS 

Having derived the magnetic field configuration, we can 
now derive the current density J from Ampere’s law. 


J = 


-Bn 


I I l± (x^u) + 1 


4irr^ I R^ 


«(« + /) dx~\ 

Bern P„(cos 0*) 


Considering the lateral balance of forces across the 
sheet, we require that 


(9) 
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where p is the gas pressure. Integration of equation (9) 
then yields an expression for the gas pressure in the 
sheet. 




8nx^ 




where Be(x) is the external magnetic field strength. This 
is the condition one would expect were there no mag- 
netic field inside the sheet. We now see that the depar- 
ture from this condition is entirely due to finite 
resistivity effects since these effects introduce an 
internal magnetic field. 

Evaluating equation (10) at x = 1 , we can solve for n 
as a function of the pressure difference across the sheet. 


where Ui =u(x - 1) and a - (£’—l)/0 •£’ = /i(l)/Pe(l)^ 
with / j(l) being the pressure in the sheet at x " 1 and 
Pe(l) the pressure just outside, and (p = IBnPe{l) 
represents the ratio of magnetic to gas pressure outside 
the sheet at the reference level. A cursory examination 
of equation (11) might suggest that 6 1 1 Rem (since 

5 1/n), yielding an extremely small value for the sheet 

thickness. However, the velocity terms Mi and 
dldx(x^u)\^ = 1 are expected to be very small. Also, if 
o - Ps(l) = i^e(l) + 5(9^ /Stt [Pneuman and Kopp, 
1971] so that a= 1. In our case a will not be exactly 
one but extremely close to one since » 1 . Hence, 
the right-hand side of equation (1 1) will be significantly 
reduced by these considerations, resulting in an 
increased sheet thickness. 

To evaluate a, Mi , d/dx{x'^u)\^ = 1 , it is necessary to 
calculate the flow along the sheet. To do this, we must 
take into account the radial jXB force, which because of 
the high current density is not small. (It will be shown 
that the radial jXB force is actually independent of the 
electrical conductivity close to the sun.) The radial 


momentum equation^ once the density is eliminated via 
the continuity equation, then becomes approximately 

( - 2 ^ 

\ u) dx X (1 + oat>y4iRgff^* 

where - GM^ffoVs^ . In deriving equation (12), we 
have assumed that 



Equation (12) without the last term is the solar wind 
equation used by Parker [1958] for an isothermal 
corona. The critical point in that case is seen to be at 
x^ = \pl2 where m ” 1 . The last term here is the radial 
component of the magnetic force, which is negative and 
hence acts to retard the expansion. As a result, the 
critical point is moved outward and occurs where the 
left-hand side of equation (12) vanishes, yielding 

2 / )p\ 


For no field, 0 = 0 and equation (13) reduces to 
Xc - 0/2 [Parker, 1958] . 

We now integrate equation (12) to obtain 


0 

—-lnu-2lnx- - 
2 X 


2n'^(p 

— - — — — -const 


which, evaluated both at the critical point and at x = 1 , 
yields 


1 

2 


0 

'2 In Xp— 


2n'^<t> 

SuiR^i^Q +tt0)x^^ 


Ml ^ 2n^0 

—--In Ml - 0 — — — — — 

2 jMii?^^(i+#) 


( 14 ) 


Equations (13) and (14) represent two coupled equa- 
tions forx^j and M l . However, the quantities n and a are 
also functions of m 1 and must be evaluated. 

To do this, consider the portion of the streamer below 
ro which leads up into the sheet. Here the flow is along 
the field and no magnetic force exists along field lines. 
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Conservation of momentum along the field then shows and, after some manipulation, equations (13) and (14) 
that become 


Psil) = Po 


[-('-?)•> 


(Mo — Ml^) 


(15) 


where is the pressure at the solar surface 

and Uq is the velocity there. Noting that Uq « , we 

can write equation (15) in the form 

Pp)=p* 

where is the pressure at the neutral point that would 
be obtained for a = «> (the hydrostatic value). Now we 
have 


^_E-1 _ 

“ <P <p 

^ [Pc/Pe(i)]exp(-MiV2) -j 

<P 

But Pc*!Pei}) - 1 + [Pneuman and Kopp, 1971] , 
Hence 



Noting that «i^/2 « |/« I we finally obtain a single 
transcendental equation forx^. 


In 




3f)X(?{2Xc-^ii) 




+ (2x0 t//) 


7 

• J: = 0 


(19) 


The numerical procedure is to solve equation (19) for 
Xc\ then Ui can be found from equation (17) or (18). 
The sheet thickness 6(-7rr/^«) can then be evaluated 
from equation (16). Figure 2 shows a plot of 6 at the 
coronal base as a function of temperature for various 
values of (j>. For reasons to be discussed later, these 
values should properly be considered lower limits to the 
actual sheet dimensions. 


(i + 0)exp(-«j V2)-/ 
a= — — ^ — — 

<t> 

and, since Ui « 1, we can expand the exponential in 
the numerator and obtain 




Tdo® “K) 


Using equation (12), we also can evaluate d/dx(x^u)\x~i 
and find 


dx 






2n'^(t> 


Equation (1 1) now reduces to 


Figure 2. Sheet thickness at the coronal base as a 
function of coronal temperature. These values are based 
upon an electrical conductivity unaffected by the mag- 
netic field and should be considered lower limits. 

We now summarize the solutions (approximate) for 
the pertinent properties inside the sheet using the 
approximations 7?^^ » 1 and 



nin + 1) /• ^ dx 

^em 
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Ba r- cos 


RemX^u 


[(-1)^%-] 


/ = 


4W(p:^ 


P= Ppix)-^ 


cos 


V 

8wc‘ 


[("4)*’* j] 

; COS=[(„4)«+i] 


small velocities in the sheet near the coronal base reflect 
(20) the retarding force of the internal magnetic field. The 
density enhancement for various field strengths is shown 
in figure 4. As expected, the enhancement close to the 
sun increases with increasing field strength (increasing 0) 
In this model, the enhancement falls to zero at large 
distances; however, this behavior may be due to the 
omission of thermal conduction effects in the model. It 
* is expected that thermal conduction will be inhibited in 
^ ^ the sheet by the transverse magnetic field, whereas con- 

duction is essentially unimpeded outside. As shown by 
Pneurmn and Kopp [1970], the inclusion of this effect 
^ 22 ) could tend to maintain significant density enhancements 
at large distances. 


The velocity u is given by the equation 
^ — Inu +2lnx^\p 


H) 


2n^<S> 


SuiRgmil +<P) 


— (^-4) 

+ <t>) \ X^/ 


(24) 


Figure 3 shows the expansion velocity both inside and 
outside the sheet as a function of radial distance. The 



Figure 3. Variation of expansion velocity with radial 
distance for (f> = fro ^2 Rq, and 1.0X10^ ""K, The 
solid curve represents the velocity inside the sheet 
whereas the dashed curve is the velocity just outside. 
Note that the flow in the sheet near the base is 
extremely small due to the inhibiting effect of the 
magnetic field. 



Figure 4. Density enhancement (ratio of density in 
the sheet to that outside) as a function of radial dis- 
tance for the cases (p = 0.1, f and 10. 


It is of special interest to note the dependence of the 
various physical quantities on the electrical conductivity 
Rem- Solutions of equation (19) over a wide range of 
magnetic Reynolds numbers reveals that the location of 
the critical point (x^) is not a very sensitive function of 
Rem- For example, for coronal temperatures in excess of 
10^ °K, Xe only varies from 1 to about 4 for 
1 < Rem <10^"^. Equation (17) then shows that 
Ui and from equation (16) we have«oc^^^. 

This means that the sheet thickness 5 varies as 
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We can now use equations (21) and (22) to evaluate the 
transverse field in the sheet (5^) and the electric current 
density J. Close to the sun,« ^Ui Hence, near 

X - 1 , we have (since n ^ ). 


Be ^ Km 


At larger distances, m is of order 1 and 


R 


Ba Qc 
O ^3 


' 2/3 

em 


x» 1 


The current / has the dependence 


R 


J oc 


em 


Consequently, the radial JXB force is indepen- 

dent if R^ffi close to the sun but varies as at large 
X, signifying that the magnetic stresses are important at 
lower levels but do not significantly affect the expansion 
at large distances. 

DISCUSSION 

Since the neutral sheet model presented in the previous 
sections is based upon the MHD approximation, it is 
appropriate now to examine its validity in structures 
whose lateral dimensions may be of the order of a mean 
free path. For example, a coronal temperature of 10^ 
and an electron density of 10^ /cm^ corresponds to a 
mean free path of 100 km, which is larger than the 
dimensions shown in figure 2. However, the Larmor 
radius, probably the more appropriate dimension in this 
case, is only 20 cm for a field strength of 1 gauss. This 
dimension is small relative to the thickness so that elec- 
trons are very effectively tied to the field lines. For this 
reason, ion-electron collisions will be frequent making 
the MHD approximation a reasonable initial approach to 
the problem. 


magnetic field. From equation (21) the transverse field 
Bq at the center of the sheet at x = 1 for Bq =1 gauss, 
ro = 2 R 0 ^ and a temperature of 10^ (/^^m ~ 

4.5X10^"^) is found to be only SXIO"”® gauss and 
decreases outward rapidly . The product of the cyclotron 
frequency and collision time for electrons in this field is 
about 70. Hence, Rem could be reduced by a factor as 
great as 2X10”^. This consideration would increase the 
sheet dimension by the factor 17 yielding widths in the 
inner corona in the order of 100 km. This figure corre- 
sponds to a thickness of 10,000 km at 1 AU. Another 
important consideration may be various forms of micro- 
turbulence and instabilities in the plasma, which would 
tend to increase the “effective” collision frequency, 
resulting in a further reduction in electrical conductivity. 

With regard to the generating mechanism for type III 
radio emission, we should note that the neutral sheet 
described here contains a weak but finite transverse mag- 
netic field. This field would have a retarding effect on 
electrons traveling outward along the sheet. One might 
well ask whether the concept of electrons moving out- 
ward at velocities of the order 0.3c across a magnetic 
field of 10“^ gauss is realistic. On the other hand, sta- 
bilization of the beam in such a configuration seems 
more possible here than in the ambient corona since the 
rapid increase in field strength outward from the center 
of the sheet should have a confining effect on the 
electrons. 
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DISCUSSION 


R, A. Kopp It’s always nice to not worry about the convergence problem by iterating 
just once, but it seems like your thickness is varying by three orders of magnitude from 
the assumption of radial. Are you sure this is conclusive? Should’! you work it out? 

G. W, Pneurmn I would say it is not a conclusive result, no. 

F. W. Perkins As far as I can see, the reason why we require finite conductivity is 
because you assume that there was flow out. I suspect there is a possibility of having a 
static solution, too, with infinite conductivity. Is the flow really observed or necessary? 

G, W. Pneurmn I believe it h necessary to satisfy the condition of almost zero pres- 
sure at infinity. 

T, G. Cowling I would like to ask, in view of the background to these latest contri- 
butions, whether in fact we are approaching the idea that the layer that you get at the 
boundary between two sectors is to be regarded as a shear layer, the equivalent of the 
earthquake boundary we’ve heard so much about? Is it a thing where the field virtually 
becomes detached, field from the one side is to be regarded as detached from the field on 
the other side? Is that a picture we have come to? And is there a genuine reason for that? 

G. W. Pneurmn I’m not sure. In this particular model I would consider it that way. I 
would like to look at it as the loops, the low down loops, actually being drawn up inside 
the sheet so that actually there are never really any open field lines inside the sheet; sort of 
an extension of the lower part of ihe^^streamer. 

Unidentified : Y6u’<f sMl^have a discontinuity involved above and below? 

G, W. Pneurmm Nov there is no discontinuity. 
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